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Abstract — In this paper, we propose a signal-selective spectrum 
sensing method for cognitive radio networlis and specifically 
targeted for receivers with multiple-antenna capability. This 
method is used for detecting the presence or absence of primary 
users based on the eigenvalues of the cyclic covariance matrix of 
received signals. In particular, the cyclic correlation significance 
test is used to detect a specific signal-of-interest by exploiting 
knowledge of its cyclic frequencies. The analytical threshold 
for achieving constant false alarm rate using this detection 
method is presented, verified through simulations, and shown 
to be independent of both the number of samples used and 
the noise variance, effectively eliminating the dependence on 
accurate noise estimation. The proposed method is also shown, 
through numerical simulations, to outperform existing multiple- 
antenna cyclostationary-based spectrum sensing algorithms un- 
der a quasi-static Rayleigh fading channel, in both spatially 
correlated and uncorrelated noise environments. The algorithm 
also has significantly lower computational complexity than these 
other approaches. 

I. Introduction 

Spectrum sensing is a key step in effectively realizing 
cognitive radio networks (CRN). In the CR access paradigm, 
secondary users (SU) in a CRN are allowed to access spectrum 
reserved for use by licensed or primary users (PU) given 
that 1) those resources are either currently unoccupied or 2) 
interference to the primary network is kept under an acceptable 
level |[T|- The main goal of spectrum sensing is to accurately 
and efficiently detect the presence or absence of a PU in a 
given band. 

Several spectrum sensing methods have been proposed in 
the literature |2 |. In general, these methods can be categorized 
as being based on either energy detection, spectral correla- 
tion (cyclostationarity), or matched filtering. Energy detection 
requires the least prior knowledge about the signal, while 
matched filtering requires the most. Spectral correlation-based 
techniques lie in between, requiring either prior knowledge 
or accurate estimation of the cyclic frequencies present in 
the PU transmission signal. Although energy detection offers 
the lowest computational complexity and is the optimal blind 
detector in the presence of i.i.d. noise, its performance relies on 
accurate knowledge of noise power due to the SNR wall phe- 
nomenon 1 3 1 . The detection performance of energy detection 
also degrades in a temporally correlated noise environment. 
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In some scenarios, spectral correlation-based methods offer 
several advantages over other spectrum sensing approaches. 
Unlike energy detection, they do not suffer from the SNR 
wall issue. These methods are also resilient to temporally 
correlated noise and enable signal-selective spectrum sensing 
where the presence a signals-of-interest (SOI) can be detected 
based on their unique cyclic features due to their modulation 
type, symbol rate, and carrier frequency f4l. 

One issue encountered with all spectrum sensing methods 
is the effect of a fading in the channel between the PU and 
SU. There is some probability of low detection performance 
whenever the channel is in a deep fade. This can be alleviated 
by exploiting spatial diversity either through the use of coop- 
erative spectrum sensing f5l or if available, the use of multiple 
antennas. As a result, spectrum sensing algorithms exploiting 
multiple antennas have received considerable interest |j6|, ||7j. 

Algorithms that leverage the cyclostationarity property have 
been applied in the past for multiple antenna receivers. In [8|, 
the sum of the individual cyclic correlation for each antenna 
was proposed. Such methods are considered post-combining 
techniques since knowledge of the channel state information 
(CSI) is not exploited. On the other hand, pre-combining 
techniques which utilize an estimate of the CSI to varying de- 
grees have been shown to have better performance. A method 
based on equal gain combining (EGC), was investigated in 
|9|. This method uses phase offset estimates to align the 
raw samples from each antenna. The aligned signals are then 
summed before finding the cyclic correlation. Finally, a Blind 
Maximal Ratio Combining (BMRC) scheme was evaluated in 
1 ,10 ] which utilized the singular value decomposition (SVD) 
to find an estimate of the CSI and applied MRC on the raw 
samples. 

In this paper we propose a spectrum sensing algorithm 
designed for use in a multiple antenna system based on the 
cyclic correlation significance test (CCST). The CCST was 
used in [11 J to perform cyclostationary source enumeration 
using an information-theoretic criterion. However, the per- 
formance of this statistic in the context of multiple-antenna 
cyclostationary spectrum sensing has not been investigated in 
prior work. The performance of this method in fading channels 
has also not been evaluated. In this paper, we show that 
the proposed method has better detection performance than 
existing methods and has less computational complexity than 
BMRC from |10|. We investigate the algorithm's performance 



in both uncorrelated and correlated noise environments. The 
scheme's robustness to intereference is also shown. 

The rest of the paper is organized as follows. The system 
model is introduced in Section [ll] including a brief discussion 
of cyclostationarity. The proposed algorithm is detailed in Sec- 
Numerical results for various scenarios are presented 
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tion 

in Section IV Finally, the paper is concluded in Section |V] 

II. Background and System Model 

A. Background on Cyclostationarity 

A signal is considered to be cyclostationary if its statistical 
properties are periodic. Equivalently, if the cyclic autocorrela- 
tion function, defined as: 



t--y-i^--^dt, (1) 



is non-zero with some r for at least one a ^ 0, the signal is 
said to exhibit second-order cyclostationary property with a 
referred to as the cyclic frequency. 

For example, in BPSK signals, cyclostationary features exist 
dX a = ijr and at a = ±2/c + where Tf, is the symbol 
period, /c is the carrier frequency, and k e { — 1,0,1}. 
Detailed analysis of the cyclostationary features for various 
digital modulations can be found in |4|. 

B. Signal Model and Assumptions 

We adopt a similar signal model as that used in pO) . 
The spectrum sensing problem is to decide between two 
hypotheses: 'Ho, where the signal is absent; and T-Li, where it is 
present. The received signal samples under the two hypothesis 
are given respectively as follows: 



x{n) 



ri{n), 

s{n) + ■q{n) 



"Ho 
Hi. 



(2) 



The received signal, sampled at a rate of l/Tg, forms 
M streams coming from each antenna with N sam- 
ples each. This received signal is defined as x(n) = 
[a;i(7T.), X2{n), . . . , XM{n)]'^ . The received signal is the super- 
position of P signal sources (including both the SOI and any 
interferers) and can be expressed in vector form as 



X {n) — ^ hj (n) (g) Sj (n) + rj {n) . 



(3) 



where (E) is the convolution operation over n and ri{n) is the re- 
ceiver noise denoted by //(n) = [771(71), 772 (?^), . . . , t/a/ (n)]^, 
where every 77^ is a purely stationary Gaussian random process 
(i?^(r) = for any a 7^ 0) with variance of cr^. For 
simplicity, we restrict that only one PU transmission, Si {n), is 
considered a SOI and that it is cyclostationary with a unique 
cyclic frequency at a ~ uq. The other sources are considered 
interferers. The channel experienced by each of the P sources 
is given by hj{n) = [hji{n), hj2{n), . . . ,hjM{n)]'^ , where 
hjk{n) is the channel between the jth source and the fcth 
antenna. We assume that the channel, although unknown to 
the receiver, stays constant over the frame of observation. 



C. Spatially Correlated Noise Environments 

In the case of spatially correlated noise, which can happen 
when there is substantial ambient noise in the band, we 
model ri{n) to have a covariance matrix given by {R^j^jlij = 
(T,,/?!*"^!. Thus with p = 0, the covariance matrix simplifies to 
an identity matrix giving spatially white noise, while p = 1 
gives fully correlated noise over all antennas. Varying degrees 
of partial correlation can be achieved by setting < p < 1. 

III. Proposed Algorithm 

In this section, we describe the proposed method. We focus 
on a single cycle frequency detection, but this approach could 
be generalized to multi-cycle detection. The key idea of this 
detection algorithm is based on the theory of canonical variates 
or theory of common factors. As discussed in | |12j , and 
subsequently utilized in pT| , the number of common factors 
between two M x 1 time-series vectors u(7i) and v(7i) can be 
inferred from the rank of the matrix 



(4) 



where the covariance matrices are defined as the time average 

R-uv = jjJ2n=o^i''^)'^^i''^)- I" non-asymptotic case, 
when this matrix is always full rank, a threshold on the 
eigenvalues can be applied to determine the asymptotic rank 
with a given confidence level. The same criterion, referred 
to as the Cyclic Correlation Significance Test (CCST), was 
applied for cyclostationary source enumeration in LllJ by 
taking u(77) = ^{n) and v(77) — x{n — T)e~^'^'""^^ for 
a given lag r and cyclic frequency a. Additionally, some 
cyclic frequencies, such as those located on a = ±2/^ for 
BPSK, only appear in the conjugate cyclic correlation matrix 
by instead taking v(7i) = x.*{n~ -^^g-jSTranT^^ 

We further adapt the CCST into the binary hypothesis test 
required for signal-selective spectrum sensing. Under T-Lo, Q 
will have zero rank as — 00. Thus, by applying a threshold 
to all the eigenvalues we can infer the presence or absence 
of the PU of interest using a M antenna receiver. Prior to 
performing the detection, we pick the lag r that provides the 
best detection performance based on the modulation format 
used by the PU. This could be done off-line by performing 
the maximization. To = argmax,- \R'^°{t)\. 

The steps of the algorithm are summarized as follows: 
1) Estimate the covariance matrix of size M x M 

^ N-l-Ta 

Rxx(to) = - J2 x(77)x^(n-To). (5) 



n=0 



2) Estimate the cyclic correlation matrix using a cyclic 
cross-correlogram at cyclic frequency a and lag tq, 
defined as 



RS°(to) 



1 



JV-l-ro 



- ^xHx>-To)e-^■2-"«"^^ (6) 

»i=0 

We will refer to the To-lag covariance matrices for both 
conventional and cyclic autocorrelation function simply 
as Rxx and R^x from this point for the sake of brevity. 



since other r are not utilized by the proposed algorithm. 
The dependence on r will be indicated explicitly when- 
ever necessary. The CCST is then calculated by finding 
the matrix 



3) Calculate its singular value decomposition (SVD) 



R^USV 



H 



(8) 



where U and V are unitary matrices whose columns 
are the left and right singular vectors, respectively and 
S = 1m [mi; M2- ■ ■ ■ ^ is the diagonal matrix of 

real, non-negative singular values, ^i, and 1m is the 
identity matrix of size M x M. 
4) Compute the test statistic by taking 



T:^ = -N\n\{{\-^Ji{). 



(9) 



The factor is used to scale the test statistic so that 
its distribution is independent of the number of samples 
used. 

5) Decision: ^tua "^^ where 7 > is a threshold chosen 
to achieve constant false alarm rate (CFAR) which will 
be discussed in the following section. 

Note that all x^(n) can be replaced with x^(n) if the 
conjugate cyclic correlation matrix is needed. We refer to these 
two metrics as the non-conjugate cyclic correlation signifi- 
cance test (NC-CCST) and the conjugate cyclic correlation 
significance test (C-CCST) respectively. 

A. Probability of False Alarm and Threshold Computation 

Two key parameters are used to evaluate the performance 
of spectrum sensing algorithms. The detection probability or 
Pd is the probability of being at Hi and accurately detecting 
the PU (Pd = Pr(7;x > 7 I ^i))- On the other hand, the 
false alarm probability, Pja, is the probability of being at Ho 
and mistakenly detecting a PU {Pfa = Pr(7^x > 7 I '^o))- 

It has been shown in fTZ| that the limiting distribution 
{N -> 00) of the test statistic based on the NC-CCST 
in ^ approaches a distribution with degree-of-freedom 
A/^. Following a similar proof, it can be shown that the 
distribution for the C-CCST is also but with degree-of- 
freedom M{M + 1). 

Based on these asymptotic distributions of under "Ho, 
the detection threshold 7 can be set to achieve a desired Pfa 
by satisfying 

/>oo 

/ f^2ix)dx^Pfa, (10) 
J ■7 



where 



\m{m 



I) 



NC-CCST 
C-CCST 



(11) 



and fxli') is the probability density function (pdf) of a x^ 
random variable with degree-of-freedom k. 
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Fig. 1 . Verification of tlie asymptotic distribution of tlie proposed test statistic 
(C-CCST) under Ho witli different number of antennas (M = {2,4}). 
Simulated pdfs for different noise variances (0-2 = {1,10}) show both the 
accuracy of the approximation and the test statistic's independence to noise 
power. The same number of samples A'^ = 4000 per antenna is used. 



These asymptotic distributions are verified to closely match 
simulation in Fig. [T]for N = 4000. Due to the scaling factor 
in (j9]l, the distribution is independent of N if the number of 
samples is large enough. The empirical pdfs for two different 
(jfj values are also shown to demonstrate how the test statistic's 
distribution under Hq is independent of noise power. 

B. Comparison With Existing Approaches 

The algorithms for multiple-antenna spectrum sensing that 
exploit cyclostationarity that are currently in the literature 
can generally be classified into two categories. Techniques 
that do not exploit any knowledge of the CSl are referred 
to as post-combining methods. The simplest method to do 
this is to find the sum of the spectral correlation test statistic 
estimated individually from each antenna ||8|. We refer to 
this approach as SUM-MSDF (where MSDF means Modified 
Spectral Density Function). The MSDF is defined as the 
spectral correlation function normalized by signal energy as 
discussed in 1 10|. 

Another existing approach is to sum the raw signals from 
each antenna and then perform a single spectral correlation 
test. However, we encounter a problem when the channel is 
not simply AWGN but instead has random fading. In this 
case, each antenna will have some unknown phase offset 
and attenuation. Thus, simply adding the raw signals non- 
coherently would decrease the probability of detection. This 
problem is addressed in |9| by first aligning the phase of 
each antenna. An estimate of the relative phase difference 
between each antenna is calculated by finding both the cyclic 
correlation of one antenna chosen as reference (auto-SCF) 
and the cross-cylic correlation of every other antenna and the 
reference antenna. The phase difference can then be extracted 



from these two. We refer to this method in our comparisons 
as EGC. 

Finally, MRC is used in pO|. Blind channel estimation is 
achieved by taking the vector corresponding to the highest 
singular value of (rob as an estimate of the channel, h. The 
raw samples from each antenna are combined using 



yin) 



h^x(n) 



(12) 



The cyclic correlation test is then performed on the com- 
bined samples y{n). This method is called MSDF with Blind 
Maximal Ratio Combining or BMRC-MSDF. It was shown to 
outperform the other techniques but at the cost of additional 
complexity due to the channel estimation and combining. One 
issue with this approach is the fact that the cyclic correlation 
is calculated twice. The first is used to blindly estimate 
the channel and the second to perform the detection on the 
combined samples. In contrast, the method proposed in this 
paper, which we refer to as eigenvalue-based cyclostationary 
spectrum sensing or EV-CSS only needs to perform the first 
part of BMRC-MSDF, the SVD, and uses the singulai" values 
themselves to infer the presence or absence of the PU. 

C. Advantages of the Proposed Algorithm 

As with other cyclostationarity-based spectrum sensing al- 
gorithms, one major advantage of the proposed method is its 
robustness to the noise uncertainty problem. Since the noise is 
assumed to be stationary and does not exhibit cyclostationarity 
at any a 7^ 0, its cyclic correlation approaches zero as 
— ?► 00. Thus, the effect of any error in the noise power 
estimate on the detection probability can be eliminated by 
taking more samples. However, in the interest of conserving 
power and arriving at a timely decision, both of which are high 
priority in the case of CR applications, we aim to minimize 
N needed to achieve a target Pd. This presents another, more 
subtle, issue related to noise uncertainty. 

In the non-asymptotic scenario, the SCF under T-Lq has been 
shown to depend on both N and the noise power cr^ fT3l. 
Therefore, the proper detection threshold is still a function of 
the noise variance. By incorrectly specifying this threshold, 
the detector could be at the wrong point in the receiver 
operating characteristic (ROC) curve. Equivalently, the target 
CFAR cannot be achieved. However, as previously discussed 
and demonstrated in Fig. [T] the proposed test statistic is 
independent of both cr^ and N . Consequently, the threshold 7 
only needs to be chosen once for a given number of antennas 
M to guarantee CFAR. This property has been shown for other 
eigenvalue-based approaches [14J . It derives from the fact that 
noise power estimation is built-in to the algorithm. 

D. A Note on Complexity 

We can make an approximate complexity comparison of 
the proposed algorithm with the best performing existing 
algorithm (BMRC-MSDF) by taking number of complex 
multiplications required for each under the same number of 
samples N . Since the cyclic covariance operation and the SVD 



are common to both algorithms they are not included in the 
analysis. 

Assuming the MSDF is calculated using an A^5-point 
Fast Fourier Transform (FFT) it requires in the order of 
A^log2(A^5) multiplications. In addition, {M + l)N multi- 
plications are needed to perform the MRC and normalization. 
Finally, the correlation in frequency uses NNs/2 multiplica- 
tions. Thus, the BMRC-MSDF approach performs in the order 
of iV(log2(A^5) + Ns/2 + M + 1) multiplications without the 
SVD and the cyclic covariance. 

In comparison, the proposed EV-CSS method finds the 
conventional covariance in addition to the same SVD and 
cyclic covariance as BMRC-MSDF, or in the order of NA'P 
multiplications. The operation R^xR-xx Q essentially 
the solution to a generalized linear system which can be 
seen as an LU decomposition requiring approximately 2M'^/3 
multiplications. Therefore, the EV-CSS approach requires in 
the order of NM"^ + 2M^ /3 multiplications in addition to the 
common operations with BMRC-MSDF. Since M is typically 
much less than both N and Ns, there is overall a significant 
decrease in complexity with the proposed algorithm. For 
example, if we take N = 4000, Ns = 128, and M = 2, 
(same parameters used in |10|), the BMRC-MSDF requires 
^296K multiplications while EV-CSS needs only ^16K mul- 
tiplications, without counting the common operations. 

IV. Numerical Results and Discussion 

In this section, simulation results are presented in order 
to compare the proposed algorithm with the various existing 



techniques discussed in Section III-B The PU is assumed to 
be transmitting a BPSK signal at a carrier frequency fc = 80 
KHz with symbol period of 25 /is. The SU is assumed to have 
M — 2 antennas unless otherwise specified. Each antenna is 
sampled at a rate fs ~ 320 kHz. For CFAR experiments we 
set Pfa = 0.1. 

In all experiments, the channel between the PU and each 
antenna of the SU, h, is modeled as a quasi-static Rayleigh 
fading channel, where the fading remains constant during the 
whole frame of iV = 4000 samples per antenna used for 
detection. This can be described using a channel vector for the 
ith frame as hj = [rie^^i,r2e-'^^, . . . , 7-j\/e-'^"^]^, where r„ is 
a Rayleigh distributed random variable of unit variance and On 
is a uniformly distributed random variable in [0, 2tt]. We ignore 
dispersive channels where each element of is modeled as a 
multi-tap filter, h.i{n), since cyclostationary spectrum sensing 
has already been shown to be robust to the effects of temporal 
correlation in prior work fTO^. 

For all algorithms, the same cyclic frequency located at 
aQ = 2fc is used. This cyclic feature has been shown to only 
be present in the conjugate cyclic autocorrelation which means 
the C-CCST statistic is used for the proposed algorithm. This 
feature is chosen because of its relative strength compared to 
other cyclic frequencies. The maximum cyclic autocorrelation 
at this cyclic frequency is observed at tq — 0. As for the 
algorithm that computes the MSDF, the frequency resolution 
is chosen to be /g/lOO. 
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Fig. 2. ROC curves comparing the proposed EV-CSS to BMRC-MSDF 
and to a post-combining method using EGC. The same number of samples 
N = 4000 and SNR = — 14 dB are used for all techniques. 
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Fig. 3. The effect of varying SNR to detection probability (P^) for different 
multiple antenna cyclostationary spectrum sensing techniques including the 
proposed EV-CSS under {Pfa = 0.1) under uncorrelated noise. The same 
number of samples N = 4000 is used. 



A. Spatially Uncorrelated Noise 

In the case where the noise between each antenna is uncorre- 
lated (Rt]t] = o',,Ij\/) and only the SOI occupies the observed 
band, the ROC curve completely determines the performance 
of each detection algorithm under a given SNR and N. A 
comparison of these ROC curves for the algorithms discussed 
are shown in Fig. |2] For these simulations the SNR is fixed at 
-14 dB, although the trends are consistent for other SNR in the 
range [—20, 0] dB where the experiments were conducted. As 
seen from this figure, the proposed method has better detection 
performance than any of the other algorithms. 

Interestingly, the method also outperforms BMRC-MSDF 



which as discussed in Section III-D has significantly higher 
computational complexity. Although this result initially ap- 
pears to be counter-intuitive, further experiments show that in 
the case of perfect CSI (not shown in figure) the MRC-MSDF 
has comparable performance to EV-CSS. Therefore, we can 
conclude that at very low SNR the blind channel estimates 
based on the SVD have large errors and the full benefit of 
MRC is not achieved. In contrast, the EV-CSS is able to fully 
take advantage of the information from all antennas because 
it does not directly use an estimate of the CSI. 

The effect of varying SNR on probability of detection is 
shown in Fig. [3] The Pfa is kept constant at 0.1. It can be 
seen that all the methods except EGC reach the P^; = 1 at 
dB SNR. Again, we see that the proposed method consistently 
has better detection performance at all SNRs compared to the 
rest of the algorithms. 

B. Spatially Correlated Noise 

The effect of varying number of samples, N, is shown in 
Fig. [4] The second set of plots also show a spatially correlated 
noise environment with p = 0.5. All methods, with the 
exception of EGC, are robust to spatial correlation and in fact. 
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Fig. 4. The effect of varying number of samples, A^, to detection probability 
(Pd) of the various techniques. Solid lines are for spatially uncorrelated noise 
environments while dashed lines are for p = 0.5. 



the proposed method shows a slight improvement in spatially 
correlated noise. We can also see that the proposed method 
shows good performance even at low number of samples while 
the MSDF based methods require more samples to achieve 
comparable performance. 

C. Effect of Interfering Signal 

We test the robustness of these algorithms in the presence 
of a strong co-channel interferer by introducing another BPSK 
signal with the same symbol rate with 30% spectral overlap. 
The effect of the interferer on the detection performance is 
shown in Fig.|5]as the signal-to-interferer ratio (SIR), defined 
as the ratio of the interferer power to SOI power, is varied 




SIR 



Fig. 5. The effect a co-channel BPSK interfer on detection probability with 
30% spectral overlap. The number of samples used for all antennas is A'^ = 
4000 and the noise level is kept constant at ct^ = 1. 

from -20 dB to dB. The noise is kept constant at cr^ — 1. 
The proposed algorithm shows very good signal selectivity, 
giving Pd close to 1 even in the presence of a co-channel 
inteferer with 100 times the SOI's power. By performing the 
correlation entirely in time domain, the proposed method is 
able to suppress the interferer much better than the MSDF. 

D. Effect of Number of Antennas 

The effect of number of antennas, M, on detection accuracy 
is studied in Fig. [6] Note that for EV-CSS, to keep the the Pfa 
constant at 0. 1 the threshold must be set to a new value based 
on ( [T0| . On the other hand, for the rest of the algorithms, the 
threshold is set for different N, M, and SNR. The SNR across 
all antennas is assumed to be the same. Since the total number 
of samples increases with more antennas {Nt = MN), we 
expect both algorithms to perform better with higher M. In 
addition, more antennas also introduce spatial diversity which 
reduces the probability of all antennas being at a deep fade 
during the sensing period. Similar to previous results, the EV- 
CSS has better performance than BMRC-MSDF for different 
values of M. 

V. Conclusion 

A multi-antenna cyclostationary-based spectrum sensing 
algorithm based on the cyclic correlation significance test 
was proposed. The method was shown to outperform current 
multiple antenna signal-selective spectrum sensing methods in 
the literature. The computational complexity of the algorithm 
was also compared with that of the best performing existing 
algorithm that uses MRC by blindly estimating the CSI 
and was shown to require substantially less multiplications. 
The detection threshold for CFAR was also determined both 
theoretically and via simulation to be independent of the 
noise variance or the number of samples. This means that 
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Fig. 6. The effect of the number of antenna on the detection probability. 
The same number of samples per antenna N = 4000 is used. 

a single threshold is required for a given number of antenna, 
eliminating the need for separate noise estimation. Future work 
includes theoretical performance analysis of the proposed 
algorithm. 
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